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ABSTRACT 

The formation of massive spheroidal galaxies is studied on a visually classified sample ex- 
tracted from the ACS/HST images of the GOODS North and South fields, covering a total 
area of 360 arcmin 2 . The sample size (910 galaxies brighter than «F775W = 24) allows us to 
explore in detail the evolution over a wide range of redshifts (0.4 < z < 1.5; median redshift 
0.68). Three key observables are considered: comoving number density, internal colour distri- 
bution; and the Kormendy relation. The comoving number density of the most massive galax- 
ies is found not to change significantly with redshift. Extrapolation of our sample to z=0 gives 
an increase of the comoving number density of M s > IO^Mq galaxies by a factor 2 between 
z=l and z=0, in contrast with a factor ~ 50 for lower mass galaxies (10 10 < M s /Mq < 10 11 ). 
One quarter of the whole sample of early-types are photometrically classified as blue galaxies. 
On a volume-limited sample out to z<0.7, the average stellar mass of the blue ellipticals is 
5 ■ 10 9 M Q compared to 4 • 10 10 M Q for red ellipticals. On a volume-limited subsample out to 
z= 1.4 probing the brightest galaxies (My < —21) we find the median redshift of blue and red 
early-types: 1.10 and 0.85, respectively. Blue early-types only amount to 4% of this sample 
(compared to 26% in the full sample). The intrinsic colour distribution correlates overall bluer 
colours with blue cores (positive radial gradients of colour), suggesting an inside-out process 
of formation. The redshift evolution of the observed colour gradients is incompatible with 
a significant variaton in stellar age within each galaxy. The slope of the Kormendy relation 
in the subsample of massive galaxies does not change over 0.4 < z < 1.4 and is compati- 
ble with z=0 values. The "zero point" of the Kormendy Relation (i.e. the surface brightness 
at a fixed half-light radius) is 1 mag fainter (in the B band) for the subsample of low-mass 
(M s < 3.5 • 10 10 M Q ) early-types. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: stellar 
content 



1 INTRODUCTION 

For several decades the formation of luminous early-type galax- 
ies has occupied a central debate in observational astrophysics. A 
comprehensive literature on these objects has convincingly estab- 
lished that the bulk of the star formation in these systems takes 
place at high redshift (z ^ 2 e.g. iBower et al.ll99 2: Stanfor d" et all 
1998) with possibly no star formation activity thereafter - a model 
that is often referred to as the 'monolithic collapse hypothesis'. 
This is borne out primarily by the optical properties of early- 
type populations and their strict obedience to simple scaling rela- 
tions over a large range in redshift. This includes the small scat- 
ter in the early-type 'Fundamental Plane' (e.g. |j0rgensen et all 
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1 19961 : ISaglia et"al 1997) and its l ack of evolution with look - 
back time (e.g. iForbes et alj|l998t Ivan Dokkum & Fraidfl996h . 
the hom ogeneity and lack of redshift evolution in their optical 
colours jBower et alj|l992l : lEllis et alj[l997l ; Istanford et ai]| 19981: 
van Dok kum et 310 2000) and evidence for short (< 1 Gyr) star 
formation timescales in these systems, deduced from the over- 
abundance of a elements(e.g. iThomas et alj 19991) . 

While the monolithic hypothesis reproduces the optical prop- 
erties of early-type galaxies remarkably well, it does not sit com- 
fortably within the currently accepted ACDM galaxy formation 
paradigm, in which early-type galaxies are thought to form through 
a hierarchical buildup process. Semi-analytical models, that com- 
bine simulations of hierarchical structure formation with simple 
recipes to describe the baryonic physics (star formation, chemical 
enrichment and feedback from the supernovae and active galactic 
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nuclei) have enjoyed broad success in reproducing several proper- 
ties of the galaxy population in the low redshift Universe (see e.g. 
ICole et al.ll200d : iHatton et alj|2003ri . 

While the predicted bulk photometric prop erties of early-type 
galaxies have met with reasonable succes s (e.g. iKavirai et alj2005l ; 
|Pe Lucia et alj200d : lKa"virai et alj2007l : lKavirail2008h . difficulties 
have remained in reproducing the short star formation timescales 
implied by the high alpha-enhancement ratios. While improved 
prescriptions may produce better agreement to the observed ratios 
(Pipino et al. 2008, in preparation) this remains one of the most 
important challenges to the current incarnation of semi-analytical 
models in terms of reproducing the properties of early-type galax- 
ies. 

Nevertheless, the discovery, using rest-frame UV data, of 
widespread low-level recent star formation (RS F) in early-type 
galaxies, both a t low redshift {Kavirai et al.ll2007l) and at i nterme- 
diate redshifts ( IFerreras & S ilk 2000; Kav irai et alj|2008h . points 
towards a more extended star formation history than previously en- 
visaged. While much of the work that underpins monolithic col- 
lapse was conducted in the nearby Universe, early-type populations 
at higher redshift are rapidly becoming available through deep opti- 
cal surveys. These surveys offer an unprecedented window into the 
star formation histories of early-type galaxies over the last 8 billion 
years, crucial for understanding their formation and evolution. 

Surveys that combine ground-based photometry with imag- 
ing from the WFPC2 and ACS cameras on board the Hubble 
Space Telescope (HST) additionally provide the angular resolution 
(FWHM ~ 0.1") to perform spatially resolved analyses of stellar 
popul ations down to phys ical projected sizes of ~ 1 kpc at z~l 
(e.g. iFerreras et al. 2005). In agreement with their low-redshift 
counterparts, recent studies that have exploited such early-type data 
at intermediate redshifts have found that the bulk of the star forma- 
tion in these objects does indeed form at high redshift. However, 
tell-tale signatures of the recent star formation found in the UV 
studies are clearly visible. A significant fraction of early-types (up 
to 20%) exhibit blue cores that are characteristic of recent star for- 
matio n jMenanteau et alj[200ll : IFerreras et al]|2005l : |Pasquali et al.1 
2006). Furthemore, spectroscopy has revealed [Oil] emission lines 
in a similar fract ion of field early-type galaxies (Schad eetal.1 19991 ; 
iTreu et alj|2002l) . 

We present in this paper an extension of the sample in 
IFerreras et all (2005) - which was restricted to CDFS - to the full 
coverage of the ACS/HST images of the GOODS North and South 
fields. However, notice that the sample presented here does not ap- 
ply the select i on ba sed on the Kormendy relation in contrast to 
IFerreras et al.|j2005l) . The only constraint applied on the sample se- 
lection is visual classification. The depth and superb angular resolu- 
tion of the images enable us to accurately determine sizes and total 
magnitudes (which translate into stellar masses). We can determine 
in a robust way the internal colour distribution and th e Kormendy 
relatio n, out to redshifts z £ 1-2. We improve on IFerreras et al] 
d2005l) by constraining the selection method to visual classification 
(i.e. removing the Kormendy Relation constraint) and by selecting 
volume-limited subsamples. 



2 THE SAMPLE 

We describe in this section the definition of our i-band selected 
(F775W), visually classified sample of early-type galaxies. The 
sample is extracted from the ACS/HST ima ges of the Great Obser - 
vatories Origins Deep Survey (GOODS; iGiavalisco et~all f2004). 



Both North and South fields are included, covering a total area of 
320 arcmin 2 equally split between the two fields. The images con- 
sist of four deep exposures through ACS passbands F435W (B), 
F606W (V), F775W (i) and F850LP (2). The images are drizzled 
with a final pixel size of 30 mas and include a weight map used 
for the estimates of photometric uncertainties. The weight maps 
also allow us to generate realistic mock galaxies used to assess the 
accuracy of the half-light radii (see appendix A). The visual clas- 
sification was performed on the version 1.0 images, whereas all 
the subsequent work: photometry, size determination and intrinsic 
colour distribution was done on the recently released version 2.0 
images (Giavalisco and the GOODS Team, in preparation). This 
version has a significant increase in the S/N of the z and the i bands. 
The methodology presented in this paper was already applied to the 
vl.O images. We find an improvement on the photometric and size 
uncertainties, but no major difference in the overall behaviour of 
the sample. 

2.1 Size and Magnitude Estimation 

We started by selecting all objects with magnitude 
MAG.AUTO < 24.0 from the GOODS ACS i F 77sw source 
catalogues (vl.l; detections based on the zfssolp band). We 
use hereafter elliptical apertures for each galaxy, adopting the 
source centre and ellipse shape (position angle and ellipticity) as 
given in the catalogues, keeping it fixed throughout the process. 
Several objects missed by the GOODS detection, presumably 
due to proximity of a bright neighbo uring source, were inse rted 
manually, using Source Extractor (Berl in & Amoutsl 1 19961) to 
determine their shape. This initial sample contains 7,462 objects 
(3,992 in GOODS-North and 3,470 in GOODS-South). 

The GOODS ACS images suffer from residual (i.e. non-zero) 
background flux, which typically reaches a level of 1/5 of the noise 
RMS in the iF775w band. This background, along with a value of 
the noise RMS, was determined individually for each galaxy as a 
single value that corresponds to the median pixel value within a 
21" x 2l"-box, applying five iterations of clipping outliers at 2.3 
standard deviations. All sources were masked in this process, us- 
ing 2.0 times the KRON_RADIUS , taken from the official GOODS 
catalogues. 

Proper masking of neighbouring gala xies is crucial f or a reli- 
able determination of the Petrosian radius (Petrosian 1976). On the 
one hand, object masks should be large enough to cover all contam- 
inating light of close neighbours. On the other hand, if neighbour 
masks are too large they could encompass a big part of the target 
galaxy, rendering unreliable estimates of the Petrosian Radius. 

After testing several ways to define the masks - e.g. using the 
KRON_RADIUS or a preliminary Petrosian radius - we adopted 
an elliptical aperture that extends out to a surface brightness corre- 
sponding to one half of the local noise RMS ("half-noise semimajor 
axis" or "half-noise SMA"). In addition, we apply an iterative out- 
lier clipping procedure in the surface brightness calculations used 
to determine the Petrosian radius, as described below. In the deter- 
mination of the half-noise SMA, initial object masks were used that 
were based on the KRON_RADIUS from the catalogues. 

For each galaxy, we define a "Petrosian semimajor axis" (here- 
after Petrosian SMA, ap c tro), i.e., in the calcul ation of the Pet- 
rosian radius, we u se ellipses instead of circles (cf. lLotz et al.l2004l ; 
iLisker et alj|20071) . As mentioned above, the elliptical shape and 

1 http://archive.stsci.edu/pub/hlsp/goods/ 
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Figure 1. Distribution of concentration, as defined in iBershadv etaD feOOOl) 
and M20 for all galaxies with ip77sw ^ 24 (grey dots). Our final - visu- 
ally classified - sample of early-type galaxies in the GOODS North and 
South fields is shown as black crosses. At the top-right corner C and M20 
correspond to a one-component distribution with a high Sersic index. To- 
wards the lower left corner, the contribution of a second component (i.e. a 
disk) makes M20 increase faster than C, resulting in a change of slope. 



source centre of each object were taken from the GOODS offi- 
cial catalogues, and kept fixed during the process. The Petrosian 
SMA was defined as the semimajor axis a at which the local in- 
tensity falls below one fifth of the average intensity within a. The 
local intensity is measured as the average intensity within an ellipti- 
cal annulus reaching from 0.9a to 1.1a, applying five iterations of 
clipping outlying pixel values at 2.3 standard deviations. Neigh- 
bouring objects were masked using the half-noise SMA, or the 
KRON_RADIUS for sources with MAG.AUTO > 24.0. For 250 
objects (3.4%), a neighbouring source was so close that its mask 
encompassed the centre of the target galaxy, rendering a Petrosian 
SMA derivation impossible. For another 72 objects (1.0%), the de- 
rived Petrosian SMA is clearly too large from visual inspection, in 
most cases caused by diffuse light or background inhomogeneities. 

To obtain the total magnitude, effective radius, and effective 
surface brightness of each galaxy, we measured the ^775 w flux 
within a — 1.5 ap ctro , and the semimajor axis containing 50% and 
90% of this flux (050 and ago, respectively). For galaxy images af- 
fected by a neighbour mask - therefore without a Petrosian SMA - 
we used the KRON_RADIUS instead of a — 1.5ap e tro, and per- 
formed no masking of neighbours. For those whose Petrosian SMA 
was found to be too large, we also adopted the KRON_RADIUS, 
but did mask neighbouring galaxies. From a comparison of cata- 
log magnitudes (MAG.AUTO) with final magnitudes (mi) for our 
early types, we find that for 3-4% of the objects, the final magni- 
tude is brighter by more than 0.5 mag than the catalog magnitude. 
Based on our preselection of galaxies with MAG_AUTO < 24.0, 
we estimate our completeness to be 95% at rrii ^ 23.5 mag. 

The flux lost due to masked areas within the aperture was 
corrected by assigning the azimuthally averaged intensity to the 
masked pixels. We then applied a correction to the estimate of half- 
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Figure 2. Comparison between our photometric redshifts and published 
spectroscopic redshifts (see text for details). Solid (open) circles correspond 
to photometric early (late) types. The comparison is restricted to galaxies 
with ip775w < 23 (with the full sample down to ip77sw < 24, the stan- 
dard deviation of the comparison is 0.13). The long (short) dashed lines 
show the sample mean (rms). The rms given in brackets correspond to the 
scatter when 3<r outliers are removed from the sample. Inset: Redshift dis- 
tribution of the complete sample (both North and South fields; thick solid 
line), using spectroscopic redshifts where available. The other histograms 
split the sample into photometric early- (dashed line) and late-types (dotted 
line). 

light radius and apparent magnitude. This correction was obtained 
from simulations of galaxies with Sersic surface brightness profiles 
and realistic noise obtained from the GOODS/ACS weight maps. 
Notice that this correction fixes two different effects: the dominant 
one in this case comes from the fact that the apparent sizes are close 
to the resolution limit of the observations. The second correction is 
required as light measured within an aperture will lose flux, and this 
loss is related to the steepness of the surface brigh tness profile (see 
e.g. iGraham & Drivedl2005l ; [Graham et all2 005). Our simulations 
show that the correction can be written as a function of two observ- 
ables: the dominant one is the estimated R e and the second order 
one is the concentration, measured as R90 /R50 (see appendix A for 
details). 

The simulations allow us to assess the accuracy of the re- 
trieved size, magnitude and apparent surface brightness. Four mock 
samples of 910 galaxies each were generated with the same values 
of «F775\v and R e as the real sample. These galaxies wer e given 
a rand om eccentricity < e < 0.6, position angle, and ISersiq 
d 19681) index 2 < n < 4. We could sucessfully retrieve the half- 
light radius with an accuracy of 9%, the ip775w-band total apparent 
magnitude within 0.05 mag, and the average iF775w-band surface 
brightness within the half-light radius to within 0.16 mag/arcsec 2 
(see appendix A for details). 



2.2 Morphological Classification 

Galaxy morphology studies (e.g. [Abraham et alj|l994l ; lLotz et all 



Ualaxy morphology studies (e.g. Abraham et al. 1994; Lotz 
120041 : iMenanteau et alj|2004l; IScarlata et al.ll2007l)~ still lack 
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liable classification method in an unsupervised way with accept- 
able levels of contamination and completeness. Currently, all these 
methods always require visual inspection. We decided to tackle the 
selection of early-types only by visual inspection. Th e Sout h field 
sample of early-types was taken from iFerreras e t al. (2005]), with 
a total of 377 early-type galaxies, out of a total of 3,470 galaxies 
at ip775w < 24 (10.9%). For the North field we follow an iden- 
tical procedure. All sources brighter than i'f775W = 24 mag are 
inspected by the four authors of this paper. We use ds!0 to view 
each galaxy in all four bands. After a first selection, we accept an 
early-type galaxy if classified as such by at least three of us. A 
second pass of visual inspection by the same people follows, this 
time restricted to the list of galaxies already labelled during the 
first pass as early-types by one or two of us. Galaxies classified 
as early-types by three of us in this second pass are added to the fi- 
nal sample. Out of 3,992 sources in the North field, we selected 533 
galaxies (13.4%), with a total in both fields of 910 galaxies. The dif- 
ferent fraction is compatible within Poisson statistics. Furthemore, 
Cosmic Variance will certainly have an effect on the discrepancy 
between the North and South fields. 

Two "non-parametric" measures of galaxy morphology are 
shown in Fig. ED fo r our sample: the conce ntration index (C 
iBershadv etaUtoOCl) . and the M 2 o parameter dLotz et alj|2004h . 
C — 5 log 7?go/i?20, where Rx is defined as the radius within 
which one has X percent of the total flux. M20 is the normalized 
second-order moment of the brightest 20% of galaxy pixels (no 
matter where they are), and thus quantifies whether these pixels are 
confined to a small region or wi dely spread. The figure, analogous 
to Fig. 1 in lFerreras et al.l |2005), shows that the majority of galax- 
ies follow a well-defined relation in the C-A/20-plane, which con- 
sists of two linear parts with different slope. The change in slope 
marks the transition between bulge-dominated and disk-dominated 
or irregular objects: for the former, the brightest 20% of pixels are 
fully located in the centre, whereas for the latter, a growing fraction 
is located outside of the central region. 



2.3 Redshift and Rest-frame Estimates 

We retrieved the latest publicly available s pectroscopic redshif t cat- 
alogues from the GOODS fie lds: FORS2 Jvanzella et al.l2008l) and 
VVDS dPopesso et al]|2008l) in CDFS and the Team Keck Red- 
shift Survev dWirth et alj|2004l) in HDF N. InlFerreras et all j2005h 
photometric redshifts were taken from Mobashe r et al. 1 20041) for 
galaxies without a spectroscopic measurement. In this paper we 
determine our own photometric redshifts both in the North and 
South fields for consistency. Our code is based on a template fit- 
ting techni que without any pr i ors an d uses a set of eight templates 
built from lBruzual & Charlotl d2003l) population synthesis models. 
Each template is a r-model, i.e. a composite stellar population 
at solar metallicity with an exponentially decaying star formation 
rate (r). We fix the formation redshift zp=3 for all templates, with 
timescales: r = {0.05, 0.5, 1.0, 2.0, 4.0, 5.0, 6.0, 8.0} Gyr. We la- 
bel these star formation histories with a n integer t={0, ■ ■ ■ , 7}. The 
templates are tranformed according to a lFitzpatrickH1999l) law for 
Galactic reddening with a colour ex cess E(B-V)=0.012 ( HDFN) 
and 0.010 (CDFS), values taken from lSchlegel et al.l d 19981) for the 
central positions of both fields. Each model is run over a range of 
redshifts out to z=3 with the available photometry estimated over a 
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Figure 3. Distribution of various observables. The histograms separate be- 
tween the North (grey) and the South (black) fields in order to illustrate the 
effect of Cosmic Variance. 



fixed aperture. The photometry comes from the GOODS/ACS im- 
ages (B,V,i,z) and is directly extracted from the publicly available 
v2.0 imag es. For HDFN we also include ground-based photome- 
try from lCapak et al.l d2004ri : which includes KPNO/4m [/-band 
and UH/2.2m HK ' data. In CDFS optical and NIR photometry is 
taken from WFI/NT T and ISAAC/VLT as u sed in the photometric 
redshift catalogue of Mobas her et alj d2004l) . 

Figure [2] shows a comparison of photometric and spectro- 
scopic redshifts for galaxies brighter than iF775w < 23. The av- 
erage (0.002) and scatter (0.091) of the discrepancy - measured as 
Az/ (1 + zs) - is comparable to the best redshift estimates avail- 
able dMobasher et"al]|2004 [20071) . If outliers at the 3a level are 
removed from the sample, the scatter reduces to 0.061. The frac- 
tion of outliers is 2.9% with no significant dependence on galaxy 
type. 

Our photometric redshift code also identifies each galaxy with 
the template that gives the best fit. This allows us to split the sample 
into "photometric" red and blue galaxies depending on whether the 
best fit corresponds to old or younger stellar populations. We com- 
pared this new classifi cation - using the tem plates presented above 
- with those done in lFerreras et al.l d2005h. which were based on 
the photometric types from Mobashe r et alj d2004l) . In this paper 
we choose the separation between photometric red and blue ellip- 
ticals at the third template (t=2): all galaxies which are best fit by 
an exponential timescale of 1 Gyr or shorter are considered "red" 
throughout this paper. With this choice, the number of "red" early- 
type galaxies in the total sample is 670 (~75%). The fraction of red 
early-types in the North and South fields is 375/533 and 295/377, 
compatible within Poisson error bars. 

In addition, we used redshifts determined for a subsample of 
our dataset, for which IF and AP have access to ACS/HST slit- 
less grism spectroscopy (G800L) from the PEARS project (Prob- 
ing Evolution And Reionization Spectroscopically; PI S. Malho- 
tra; HST Proposal ID # 10530) The low-resolution spectra (R~ 
50 — 80) allowed us to secure redshifts with comparable accu- 
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Figure 4. Redshift evolution of the comoving number density: The B-band absolute magnitude (left) and stellar mass (right) is compared with V (< z), i.e. the 
comoving volume within a given redshift (the redshift appears in the top axis). Vertical lines delimit redshift bins used in the analysis of the number density in 
the next figure, with the upper and lower limits shown as shaded regions. Galaxies are shown as grey dots (photometric early-type) or open circles (late-type). 
The curved, thin solid lines in the left panel are the limiting «F77sw = 24 for a photometric early- or late-type. The thick line and dashed lines above and 
below illustrate the sample se paration into "brig ht" and "faint" subsets (or high and low stellar mas s in the right panel; s ee text for details). Grey areas in the 
left panel are L* values from[llbert et al. 12006). The dashed line in the mass plot follows M, from lFontana et all j2006t) . 



racy to the spectroscopic ones because of the prominent 4000A 
break found in most of these galaxies. For details we refer the inter- 
ested reader to the paper on the analysis of the grism data from the 
GOODS early-types (Ferreras et al., in pr eparation), or to a similar 
study on the HUDF dPasquali et aT]|2006l) . Out of the total sample 
of 910 early-type galaxies in this paper, we have spectroscopic or 
grism-based redshifts for 541 galaxies. Hence, 59.5% of the total 
sample have spectroscopically accurate redshifts. 

The redshifts are used to determine rest-frame properties of 
the galaxies. The physical projected half-light radius is determined 
by using the angular diameter distance for a concordance cosmol- 
ogy, namely Q m = 0.3, SIa = 0.7 and Ho — 72 km/s/Mpc. The 
best fit template is compared to the apparent total magnitude to es- 
timate the absolute magnitude and stellar mass - a Chabrier(2003) 
Initial Mass Function is used throughout this paper. From the ACS 
B (where available), V, i, z photometry we choose the passband 
with the lowest K-correction for the transformation into absolute 
magnitude. As an estimate of the amount of correction applied, 
the RMS of the K-corrections for the full sample is 0.42 mag for 
Ms and 0.27 mag for My. Regarding stellar mass estimates, even 
though the age and metallicity distribution of a composite stellar 
population cannot be constrained with broadband photometry, the 
stellar mass can be reliably constrained to within 0.2 — 0.3 dex 
provided the adopted IMF gives a n accurate representation of th e 
true initial mass function (see e.g. iFerreras, Saha & Burlesll2008h . 
In order to assess the effect of the passband on the determination 
of stellar mass we compare our stellar masses - based on iF775w - 
with the K a photometry of the CDFS field from ISAAC/VLT (Ret- 
zlaff et al., in preparation). Over the redshift range probed by our 
sample, i-band and NIR mass estimates differ by less than 0.15 dex 
(see appendix B). 

Figure[3]shows the distribution of stellar mass, projected phys- 



Table 1. The sample of early-type galaxies («F775W < 24) 



Sample 


Number 


fa 
J B 


z b 




Mi 

IVly 


R^(kpc) 


HDFN 


533 


0.30 


0.78 


4.2 ■ 10 10 


-20.9 


1.10 


CDFS 


377 


0.22 


0.64 


3.3 ■ 10 10 


-20.6 


1.26 


All 


910 


0.26 


0.68 


3.7 ■ 10 10 


-20.8 


1.15 


Sl c 


411 


0.04 


0.85 


1.2 ■ 10 11 


-21.9 


2.11 


S2 C 


280 


0.11 


0.56 


3.4 ■ 10 10 


-20.6 


1.55 



a Blue early-type fraction: Red/Blue types selected from photometric red- 
shift best type (Red: t^ 2, i.e. star formation timescale r ^ 1 Gyr, see text 
for details). b Median value. c SI and S2 are volume-limited subsamples, 
defined in §4. 



ical size, absolute magnitude and redshift. The sample is shown 
separately for the North (grey) and South (black) fields. Both fields 
have a similar range in mass or size, but the redshift distribu- 
tion illustrates the effect of cosmic variance over the individual 
160 arcmin 2 areas. The variation in the total number of early-type 
galaxies (533 in the North vs. 377 in the South) is also caused by 
such effect (the selection method being identical in both fields). 
Some of the basic properties of this sample is shown in table [T] 
The table also shows the properties of a couple of volume-limited 
subsamples (SI and S2) defined in §4. 



3 NUMBER DENSITY EVOLUTION: LUMINOSITY VS. 
MASS 

A simple way to assess the redshift evolution in the number den- 
sity of early-type galaxies is illustrated in figure [4] where the B- 
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Figure 5. Redshift evolution of the top end of the mass and luminosity 
distribtution as illustrated in the previous figure. The black dots correspond 
to the brighter (top panel) or more massive (bottom) subsample. The fainter 
subsample is shown as grey crosses. Error bars include Poisson noise and 
the uncertainty in the K-correction. The lines give fits to a power law, with 
index a, which is labelled for all subsamples. 



band absolute magnitude (Jeft) or stellar mass (right) is presented 
with respect to the comoving volume enclosed within the redshift 
of each galaxy. Hence, counting galaxies within cells in this figure 
gives the comoving number density at a given redshift. The thin 
solid lines track the flux limit imposed by our search (jab < 24). 
This limit depends on the K-correction adopted, thereby we show 
in the figure the absolute magnitude limit taking the extrema of the 
templates used in the photometric redshift analysis (i.e. types t=0 
and t=7). Galaxies are shown according to the convention adopted 
throughout this paper, namely red galaxies as grey solid dots (pho- 
tometric type t^ 2) and blue galaxies as black open circles (t>2). 
In both cases - absolute magnitude and stellar mass - we split the 
sample in two bins. In absolute magnitude we define 1 mag bins 
(dashed lines) above and below a typical L* galaxy represented 
by the thick solid line, which corresponds to M% ~ —21. In or- 
der to track its evolution with redshift, we assume an exponen- 
tial star formation history with formation redshift zp=5, timescale 
tsf = 0.5 Gyr and solar metallicity. For comparison, the shaded 
boxes delimit the value of the M% Schech ter fits of the VI MOS- 
VLT Deep Survey for early-type galaxies (Ilbert et al. 200f|). We 
decided to choose this binning criterion instead of a fixed absolute 
magnitude in order to overcome the passive fading of the stellar 
populations. The shaded areas show that our low-luminosity bin 
roughly maps L» galaxies. 

Regarding stellar mass, we choose a fixed reference mass 
(10 11 Mq; thick solid line). The sample is split in 1 dex bins above 
and below this reference (dashed lines). For referen ce, the dotted 
line tr acks the characteristic stellar mass scale of iFontana et al.l 
d2006h . A typical error bar is shown in both panels. The vertical 
lines in both panels mark the redshift bins, with the full range for 
a complete volume-limited sample given by the hatched vertical 
regions. 



Figure [5] shows the evolution in comoving number density 
with respect to absolute magnitude (top) or stellar mass (bottom) 
obtained from the previous diagram. The error bars are estimated 
from a sum in quadrature of the Poisson noise in each bin and the 
uncertainty in the absolute magnitude or mass. This uncertainty is 
computed using a Monte Carlo algorithm in which noise is added 
to the absolute magnitude or stellar mass of each galaxy. 5,000 real- 
izations are generated and each is analysed in the same way as our 
original sample. The error bars give the standard deviation within 
the ensemble. The evolution of the number density can be mod- 
elled by a simple power law n cx (1 + z) a . In luminosity there 
is no significant difference between the "bright" and "faint" sub- 
samples. However, the most massive galaxies show a remarkable 
weak trend with redshift, with a slope significantly different to less 
massive galaxies. With the adopted fits to the data, an extrapolation 
to lower redshift yields a decrease of a factor ~ 50 in the number 
density of early-type galaxies between z=0 and z=l for the mass 
range 10 10 — 1O 11 M0; whereas the decrease is just a factor 2 for 
galaxies with stellar masses between 1 11 and 10 12 Mp,. Th is result 
is presented in more detail in a letter l lFerreras et al.l 2008) where 
we show tha t this result is in agree ment with recent work in other 
surveys (e.g. IConselice et al . 2007) a nd with recent se mi-analytic 
models of galaxy formation (e.g. lKhochfar & SillfeOOot ). 



4 INTRINSIC COLOUR DISTRIBUTION 

The internal colour distribution of a galaxy can reveal important in- 
formation regarding its star formation and assembly history. Early- 
type galaxies have smooth surface brightness distributions, which 
make them ideal targets for this sort of study in GOODS. For each 
colour map three pieces of information can be gathered in a robust 
way: the average colour within an aperture (the projected half-light 
radius in our case); the radial gradient of the colour distribution, 
and the scatter about the best fit (assumed to be linear between 
colour and log R). W e follow the same technique as described in 
iFerreras et al.l J2005l) . briefly outlined here. 

The colour map is generated by convolving the image taken 
in one filter with the Point Spread Function (PSF) of the other fil- 
ter used to define the colour. This is the most robust way to reduce 
the colour trends caused by a different resolution between images 
in different bandpasses. We used a PSF gene rated from stellar im- 
ages in the same dataset. lFerreras et al. (2005) includes an appendix 
that shows the - negliglible - difference in the colour gradients and 
scatter when using either a stellar i mage from t he same data or a 
synthetic PSF generated by TinvTir rlKristl Jl993l) . 

The images span a wide range of signal-to-noise ratios on a 
pixel-by-pixel basis. Hence, it is best to bin the pixels although pre- 
serving spatial information as much as possible. A Voronoi tessella- 
tion method was adopted, in which a target S/N is specified so that 
the algorithm merges nearby pixels into tiles, keeping these tiles as 
round as possible, following the algorithm bv lCappellari & Copinl 
( 2003 ). The proc e ss is fi nished when a target S/N per tile is reached. 
In IFerreras et al.l d2005l) we showed that estimates of radial colour 
profiles are much improved as the outer parts of the galaxy have a 
significantly lower noise after a Voronoi tessellation. We impose 
a target S/N=10 per tile for each colour. The radial distribution 
of colour vs log R/R e is fit to a straight line. We also determine 
the colour scatter about this best fit by using a biweight estimate 
jBeers et al]|l990l) . 

The left panel of figure [6] shows the observer frame V — i 
colour distribution as a function of redshift, half-light radius and 
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Figure 6. Left: The observer-frame colour gradient (top), scatter (middle), and integrated colour within the half-light radius (bottom) are shown with respect 
to three observables: redshift (left); apparent half-light radius (middle) and 2F775W apparent magnitude (right). The full sample is shown in grey. Black points 
correspond to average values of subsamples (after a 5cr clipping) binned at constant number of galaxies per bin, and separated with respect to the best fit 
populations: old (filled dots; photometric type t^J 2) and young (hollow dots; t>2). Right: Redshift evolution of the intrinsic distribution of the three available 
colours, as labelled. In both figures the lines in the bottom panels correspond to the evolution of a set of exponentially decaying star formation histories with 
solar metallicity, started at zp=5 with a decay timescale of 0.5 (solid line) 1 (dashed line) and 8 Gyr (dotted line). 



apparent i band magnitude of the whole sample (grey dots). The 
larger dots correspond to binning the sample at a fixed number of 
galaxies per bin, separating them in photometric type between red 
(t^ 2; solid circles) and blue (t>2; hollow dots). Error bars rep- 
resent the standard deviation within each bin. For comparison, the 
lines at the bottom-left corner track exponentially decaying star for- 
mation histories at solar metallicity with formation redshift zf=5 
and timescales r = 0.5 Gyr (solid), 1 Gyr (dashed) and 8 Gyr 
(dotted). 

In the top panels, the radial colour gradient is shown. Galax- 
ies with red (blue) cores have negative (positive) gradients . We ob- 
tain a consistent result with respect to lFerreras et alj J2005I) . Galax- 
ies with blue cores are almost always photometrically classified as 
"blue" and galaxies with red cores are classified as "red". This is 
not a trivial statement, as the overall colour of the galaxy need 
not correlate with the radial gradient. The scatter does not corre- 
late with redshift, but there is a clear trend with size or apparent 
magnitude. This can be caused by a combination of physical and 
pas sband- shifting effects and we defer the discussion to a later fig- 
ure. 

On the right side of figure[6] a similar plot is presented for the 
three available colours from the ACS images. The B — i,V — i and 
i — z colour distributions are shown with respect to redshift. The 
lines in the bottom panels correspond to the same models as in the 
left side of the figure. The trend is consistent with all three colours, 
with blue galaxies having positive colour gradients. 

In order to extract meaningful trends of the intrinsic colour 
distributions we have to define volume-limited subsamples that will 
remove biases inherent to flux-limited samples. FigureUJillustrates 
our selection process. The complete sample is shown in absolute 
magnitude vs redshift space. The distribution shows the character- 
istic boundary caused by the iF77sw = 24 limit. The solid line 



tracks a iF77sw = 24 galaxy as a function of redshift, assuming a 
typical star formation history for an early-type galaxy (photometric 
type t=l). A compromise has to be made when generating volume- 
limited samples: either we probe as deep as possible in redshift - 
cutting in absolute magnitude, and thereby sacrificing faint galax- 
ies, or we probe as faint as possible the luminosity function - cut- 
ting in redshift. In order to explore these two choices we select two 
subsamples, shown as black dots in each panel: SI (412 galaxies; 
Jeft) is a sample targeted to explore the redshift evolution, and is 
restricted to galaxies brighter than Mv < —21. We also apply a 
cut in redshift (z<1.4). Subsample S2 (280 galaxies; right) targets 
a wide range of luminosities out to z<0.7. In this case the volume- 
limited sample extends to fainter galaxies Mv < —19. The bottom 
panels of figure UJ show the histograms in stellar mass (Jeft) and 
projected physical size of each subsample. 

Figure[8]focuses on the redshift evolution of the colour distri- 
bution and shows the observer-frame V— i colour gradient (top) and 
scatter of subsample SI. For comparison, the colo ur gradients ob- 
served in a sample of moderate redshift clusters bv lTamura & Ohtal 
( 2000) is shown as stars. Our sample is shown with the usual coding 
with respect to photome tric type: red (grey do ts) and blue (open cir- 
cles). In the appendix of lFerreras et alj d2005l) simulations of galax- 
ies with the same characteristics as the GOODS galaxies were done 
to show that from photometric errors one could expect a typical 
colour scatter a(V — i) 0.03 mag, with a maximu m of 0.05 mag 
for the faintest galaxies. Figure 15 in the appendix of lFerreras et al.l 
(2005) shows that both the observed scatter and slope are intrinsic. 

We include in the figure a simple model prediction which uses 
the local B — R colour gr adients of ellip tical galaxies in the lo- 
cal Universe as constraint JPeletierJ 1 1 99Ct) . The set of two nearly 
horizontal lines - just below the A( V — i) /A log(7?) = axis - 
assumes the z=0 intrinsic colour gradient to be caused by a pure 
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Figure 7. Volume limited subsamples of early-type galaxies in GOODS 
N+S fields. SI focuses on massive systems out to redshift z<1.4, whereas 
S2 limits the redshift out to z<0.7 in order to maximise the range of 
masses/luminosities. The top panels show V-band absolute magnitude vs 
redshift, with the iF775\y < 24 flux limit shown as a solid line. The bottom 
panels show the histograms of each subsample in stellar mass and projected 
size. 
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Figure 8. Intrinsic V — i colour distribution of subsample S 1 : The slope 
(top) and scatter of the radial colour gradient is shown with respect to red- 
shift. The sample is separated into photometric early- (grey dots) and late- 
types (open circles) . Stars are gradient me asurements in a sample of clus- 
ter early-types from Tamura & Ohta 12000). Two alternative models for the 
evolution of the colour gradient are shown, depending on whether the gra- 
dient is caused by metallicity (thin, nearly horizontal lines) or age (thick, 
curved lines). Two formation redshifts are considered for each case: zp=5 
(solid lines) or 2 (dashed lines). The dotted line corresponds to a metallicity 
sequence with formation redshift zp=5, adding to the core 10% (in mass) 
of a stellar population formed at z=0.8. 



metallicity sequence, i.e. we fix an age throughout the galaxy and 
compute the difference in metallicity that will give the observed 
colour gradient in local early-type galaxies. Then we evolve this 
model backwards in time to assess the evolution in the colour gra- 
dient. Standard models with exponentially decaying star formation 
histories are used for the fit. The timescale is fixed at r = 1 Gyr and 
the formation epoch is either zf=2 (solid line) or 5 (dashed line). 
The second sets of models - which curve down at z<0.6 - assume 
a pure age sequence for the colour gradient at fixed (solar) metallic- 
ity, with a similar set of models (in this case the bluer outer region 
is modelled by a slightly longer star formation timescale). As ex- 
pected, this model gets very steep with increasing redshift because 
of the strong dependence of colour with age. A third model (dot- 
ted line) overlays on top of a metallicity sequence formed at zF=5 
a second, younger stellar population in the core, formed at z=0.8. 
This younger component contributes only 10% in mass, illustrat- 
ing the effect that a small episode caused by a minor merger would 
have on the colour gradient. 

This result illustrates that blue cores have a radial range of 
ages, being youngest at their centres. The large blue colour gradi- 
ents observed in many of the blue cores cannot be explained by an 
"inverted" gradient in metallicity, i.e. more metal-rich populations 
outwards. If we were to explain the large range in colour between 
the inner and the outer parts of the galaxy by metallicity alone, one 
would have to invoke unrealisticaly low metallicities at the galaxy 
centre. Furthermore, this scenario would predict no redshift depen- 
dence on the fractional contribution of blue cores, in contrast to the 
observations. In this subsample (SI) of luminous {My < —21) 
early-type galaxies, the median redshift of blue and red early-types 
is 1.10 and 0.85, respectively, i.e. bluer early-type galaxies (at least 
at the massive end) contribute more significantly at z ^ 1 . 



The relation between colour distribution and galaxy proper- 
ties is explored in figure |9j where the colour gradient (top) and 
scatter is shown with respect to stellar mass; projected half-light ra- 
dius and V-band absolute magnitude for subsample S2. Blue cores 
typically have stellar masses below 1O 1O M0 and half-light radii 
R e .$ 1 kpc, a signature of downsizing. In this subsample the av- 
erage stellar mass of the blue ellipticals is 5 • 1O 9 M0 compared to 
4 ■ 10 10 Mq for red ellipticals. This strong mass-dependence sug- 
gests that AGN activity alo ne is not enough to explain the blue 
cores (Menanteau et al. 20o3)- 



5 THE KORMENDY RELATION 

The Kormendy relation (hereafter KR; Kormendvl l 19771) is a tight 
linear correlation between average surface brightness - typically 
measured within the half-light radius (R e ) - and the logarithm 
of R e . It is one of the p rojections of the Fundamental Plane 
(Diorg ovski & Davisll 19871) . which relates surface brightness, size 
and velocity dispersion. The Fundamental Plane and its projections 
are a consequence of the dynamics of early-type galaxies and their 
underlying stellar populations. KR is therefore the "next best thing" 
to the Fundamental Plane when measurements of velocity disper- 
sion are not available. The observed dispersi on of the KR at z=0 
is ~ 0.4 mag arcsec~ 2 jLa Barbera et al.ll2003h . The redshift evolu- 
tion of the Kormendy relation can be used to constrain the assembly 
and formation history. A significant change in slope with redshift 
can reveal the presence of star formation or the rearrangement of 
the stars because of "dynamical activity" (i.e. mergers). 

Figure [10] shows the Kormendy relation of both volume- 
limited samples SI and S2 (grey dots). The surface brightness is 
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Figure 9. Intrinsic (V — i) colour distribution of subsample S2: The slope 
(top) and scatter of the observed V — i colour distribution - roughly rest- 
frame B — V over the redshift range of S2 - is shown with respect to 
stellar mass (left); projected half-light radius (middle) and \/-band absolute 
magnitude (right). The black points with error bars correspond to average 
values of subsamples separated with respect to photometric redshift type: 
old (filled dots; type t^ 2) and young (hollow dots; t>2). la error bars 
are shown (after a 5a clipping). The binning is done at constant number of 
galaxies per bin. 



computed as the average within the half-light radius. Open black 
dots correspond to the upper quartile in redshift (SI, left) or stel- 
lar mass (S2, right), as labelled. The solid line is the loca l re- 
lation observed in the Coma cluster dj0rgensen et alj[l995l) and 
the crosses are a sample of 1 7 galaxies from the K20 survey 
ddi Serego Alighieri et al.l [2005h . Sample SI - which extends to 
higher redshift than S2 - has a significantly brighter zero point, 
caused by the effect of lookback time on the stellar populations. 
For comparison, a simple stellar population formed at zp=5 and 
solar metallicity undergoes a fading of ~ 1 mag in the rest-frame 
B band between redshifts z=l and 0, similar to the offset found 
in the left panel of figure [10] Sample S2 (right) has a KR com- 
patible with local early-type galaxies (again with a significant shift 
of the zero point caused by lookback time). As expected, galaxies 
with a higher stellar mass content appear brighter than less massive 
galaxies with the same size, roughly a direct c onsequence of a non- 
homo logous surface brightness distribution ( Graha m" & Guzmanl 
2003). An estimate of the typical error bar is shown in each panel. 
For the surface brightness we use the scatter of the K-corrections 
applied to obtain B-band magnitudes (0.4 mag). The true error 
should be most certainly smaller than this K-correction. The er- 
ror in the size is taken from the worst case scenario when retriving 
half-li ght radii from simulations (see appendix). The K20 sample 
from ldi Sere go Alighi eri et alj J2005) (crosses) span a narrow red- 
shift range (0.9 • ■ • 1.3), and are compatible with our high-redshift 
(Jeff) or massive (right) subsamples. 

Notice t he major difference b etween the selection criterion of 
this work and lFerreras et al.l d2005h . In that paper, a Kormendy se- 
lection was imposed on the visually classified sample, so that early- 
types that were found not to evolve into the local KR were rejected 
from the sample. In that case, a large number of low-redshift, blue 
early-types were rejected. It was discussed whether misclassifica- 
tion of bright knots of star formation or an active nucleus could ex- 
plain this rejection. In this paper we use visual classification as the 
only selection criterion. We later define volume-limited subsam- 
ples. In this process we end up rejecting similar type of galaxies as 
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Figure 10. Kormendy relation of samples S 1 (left) and S2 (right) in the rest- 
frame B-band (grey dots). Open circles dots are the upper quartile of the 
redshift (SI) or the stellar mass distribution (S2), as labelled. A typical error 
bar is shown. The data have been corrected for the (1 + z) 4 cosmological 
dimming. The line represents the local relation observed in the Coma cluster 
I Jorsensen et al. 1995). The 17 galaxies from the K20 sample are shown as 
crosses i di Serego Alighieri et al. 2005). 



in lFerreras et alj d2005l) for different reasons. Figure [4] shows that 
most of the faint, low-redshift galaxies in our sample are removed 
from SI or S2 to enforce completeness. Those are mainly blue 
early-types with small sizes. The addition of those galaxies to SI or 
S2 would result in significant scatter towards the bottom-left cor- 
ner of the Kormendy relation in figure[l0] in a similar fashion to the 
presence of the "dwarf galaxy branch" presented in the Kormendy 
Relation plot of ICapaccioli et al. I dl992h . Hence we conclude that a 
significant fraction of the galaxies rejected in lFerreras et alj ( 2005) 
because of the KR constraint must be dwarf ellipticals at z ^ 0.4. 

Figure [TT] shows the redshift evolution of the KR of subsam- 
ples SI and S2 in terms of slope (top); average surface brightness 
(middle) and scatter (bottom; using a biweight estimator). Sample 
S2 is split at the median value of stellar mass (3.5 • 1O 1O M ). The 
fit is done with a standard least squares algorithm, binning the sam- 
ples at constant number of galaxies per bin (hence the unequal sep- 
aration between dots in the figure). The arrow in the top panel rep- 
resen ts the z=0 slope of the Kormendy Relation (La Barbe raet al.1 
[2003(). Sample S 1 (black dots) - focussing on bright galaxies out to 
the highest redshifts - shows a remarkable agreement with the lo- 
cal ob servations, and with moderate redshift clusters |Ziegler et"al] 
1999). We extend the conclusion o f lZiegleretal] fl999) of a neg- 
liglible change in the slope of the KR to a wide range of redshifts 
and to a field environment. The most massive galaxies in sample 
S2 (grey triangles) also agree with the local observations, whereas 
a selection of the galaxies with the lower masses in S2 (grey stars) 
have larger scatter and steeper slopes. 

In the middle panel we compare the observed trend of the aver- 
age surface brightness with the fading of a typical old population as 
labelled in the figure caption. Notice that the least massive galaxies 
of sample S2 are consistently brighter than this estimate, reflecting 
the presence of younger stellar populations. It is interesting to no- 
tice that the ~ 0.5 mag arcsec -2 scatter about the best fit (bottom) 
does not change significantly with redshift or sample. 

To summarise, we find the Kormendy relation of z ~ 0.3—1.2 
to be the same as at z ~ 0. The slope and scatter of the correlation 
does not change and the zero point, i.e. the typical surface bright- 
ness changes in agreement with pure passive evolution. Only lower 
mass spheroidals (M s < 3. 5-10 10 Mq) show a significant departure 
from the average z — Kormendy relation, appearing brighter in 
{Hb}° and with a steeper slope. We believe this result is caused by a 
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Figure 11. Redshift evolution of the Kormendy relation. Sample SI is 
shown as black circles, sample S2 (grey) has been split with respect to the 
median stellar mass within that sample (3.5- 10 10 Mq). The bins in redshift 
are chosen such that the number of galaxies per bin is constant for each set. 
A Kormendy relation is fit, and this plot gives the slope (top), zero point 
(middle; in magarcsec~ 2 ) and biweight scatter about the best fit (bottom; 
in magarcsec ). The arr ow in the top panel is the z~0 slope obtained 
bvlLa Ba rbera et all d20O3l) and the crosses are the cluster observations of 
Ziegle retalJh999h . The line in the middle panel corresponds to the evolu- 
tion of an exponentially decaying star formation rate started at zp=5, with 
timescale r = 0.5 Gyr and solar metallicity. 



combination of younger stellar populations and a significantly dif- 
ferent surface brightness profile compa red to more massive early - 

case dGraham & Guzmarjiooa) . 



type galaxies, similar to the z 



6 SUMMARY AND CONCLUSIONS 

We present a sample of visually-classified early-type galaxies over 
the 360 arcmin 2 area covered by the AC S/HST im ages of the Great 
Obse rvatories Origins Deep Survey (GOODS, iGiavalisco et al.l 
120041) . including both the North and South fields. The classifica- 
tion is only limited in flux (iab < 24), i.e. no colour cuts or 
surface brightness profile constraints are imposed in the selection. 
The sample comprises 910 early-type galaxies, for which we com- 
pile spectroscopic quality redshifts from publicly available datasets 
for 541 galaxies (i.e. ~ 60% of the total sample). For the rest 
we estimate photometric redshifts from available data (both HST 
and ground-based). We find the comoving number density to de- 
pend strongly on stellar mass, such that the number of galaxies 
with M s > 10 U M Q only increases a factor 2 between z=l and 
z=0, whereas the number of galaxies with stellar mass 10 10 < 
M s /Mq < 10 11 increase a factor ~ 50 in the same redshift range. 
In the full, flux-limited sample, we find 26% of our early-types are 
blue (i.e. corresponding to a star formation history with timescale 
longer that 1 Gyr, or with a significant amount of young stars). This 
fraction is similar to the number of low-red shift early-types wi th a 
significant amount of recent star formation dKavirai et alj|2007r) . 
The depth and superb resolution of the ACS images enable us 



to explore the intrinsic colour distribution. We perfo rm a technique 
based on a Voronoi tessellation (already applied in iFerreras et al.l 
1 20051) that allows a robust determination of the radial colour gra- 
dient and scatter. In figure [6] we show that the colour gradients are 
strongly dependent on the overall photometric type, such that red 
galaxies have mostly red cores and blue early-type galaxies have 
blue cores. This is a non-trivial statement, as one could imagine star 
formation in the outer parts of the galaxy, resulting in blue galax- 
ies with red cores. We find central star formation in otherwise red 
galaxies is very rare. On the other hand, the majority of galaxies 
with (recent) star formation exhibit most of this star formation in 
their centers. These observations suggest an inside-out formation 
process in these galaxies, ruling out recent attempts at explaining 
elliptical galaxy f ormation with an outside-in formation scenario 
(Pipino etal]|2006|) . An inside-out scenario implies that regardless 
of the way star formation is triggered - primordial collapse ver- 
sus merging/interactions, gas is accumulated towards the centre, 
starting star formation in the central regions, subsequently "diffus- 
ing" outwards. One could argue whether the me chanism explain- 
ing these blue cores reduces to minor mergers (Elli s et al]|200ll : 
iKavirai et ai1 l2007b). Robust confirmation of the actual cause of 
these blue co res will require the an alysis of abundance ratio gradi- 
ents (see e.g. IFerreras & Silkl2002h . 

Figure[8]shows that the loca l observations of radial colour gra- 
dients in elliptical galaxies (e.g. Peletier 1990) mostly correspond 
to a range of metallicities. This can be throug ht of as the intrinsic 
equiva lent of the pure metallicity sequence of iKodama & Arimotol 
dl997l) , used in that paper to explain the colour-magnitude rela- 
tion among galaxies. Our results are in agreement with the sample 
of cluster galaxies at low an d moderate redshift dTamura & Ohtal 
l200rA lLa Barbe ra et al .120041) , but our sample increases by an order 
of mag nitude the number and a ccuracy of the colour gradient es- 
timates. La Bar bera et al.l (12004) only find optical colour gradients 
compatible with zero slope. We do find non-zero, albeit small, opti- 
cal gradients and a significant trend with stellar mass, such that the 
most massive galaxies have almost always red cores. 

Defining volume-limited samples, we find that blue galaxies 
only contribute 4% to the total number of bright (My < —21) 
early-type galaxies out to z< 1.4. Most of the blue early-type galax- 
ies in our - visually classified - sample correspond to low-mass 
systems, and are removed when taking volume-limited samples. At 
face value, we could state that significant star formation in early- 
type galaxies only takes place for galaxies less m assive than M a ~ 
10 10 M Q even out to redshifts z~l. As shown bv lBell et all d2007h . 
star formation in a substantial fraction of blue cloud galaxies needs 
to be shut off between z=l and z=0, making them evolve to the 
red sequence, thereby explaining the significant growth in mass 
there. However, in addition to considering the contribution with 
respect to mass and star formation history, it is also well known 
that blue cloud and red sequence differ morphologically, with spi- 
ral and irregular galaxies dominating the blue cloud whereas red se- 
quence galaxies are mostly early-types. Since our sample of early- 
types was selected on pure morphological grounds, th ose early- 
types w ith (recent) star formation confirm the scenario o f lBell et al.l 
d2007l) , in that these objects will not only evolve onto the red se- 
quence in terms of colour, but also in terms of morphology. 

The Kormendy relation was explored in two volume-limited 
samples and a remarkable consistency is found with the KR at z=0, 
with no significant change of the slope of the Kormendy relation out 
to z~1.2. The scatter about the best fit does not change significantly 
either. Only the average surface brightness is found to change in the 
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Kormendy relation, with an evolution consistent with pure passive 
evolution of old stellar populations. 
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APPENDIX A: EXTRACTION OF SIZES AND MAGNITUDES FROM THE ACS IMAGES 

As discussed in §§2. 1 it is important to correct for the loss of light expected for a typical surface brightness distribution. In order to estimate 
the accuracy of the retrieved size, apparent magnitude and surface brightness we ran a number of simulations, generating postage stamps 
with the same noise c haracteristics as the ones taken from the ACS images. The simulations correspond to smooth distributions of galaxies 
with a Sersic profile < jSersidll968h . The simulations realistically map the error in each pixel using the weight maps of the ACS images to 
determine a correlation between surface brightness per pixel and S/N. 

For a given choice of apparent magnitude and half-light radius we generate 501 x 501 stamps with the same pixel size as the ACS images 
(0.03 arcsec). For each galaxy in the ensemble we take a random value for the Sersic index (between n = 2 and 4); ellipticity (between 
e — and 0.6) and position angle. A total of 16 ensembles are generated for four different values of apparent magnitude (ip775w = 
{18, 20, 22, 24}) and half-light radius (R e = {5, 8, 10, 20} pixels). Each model galaxy is convolved by the PSF of ACS in the i band 
(obtained from a set of stellar images as discussed in §4) and is passed through the same code as the original sample to determine the 
half-light radius, magnitude and surface brightness. 

Figure[AT]shows the result (la error bars). Notice the systematic trend in the offset with respect to half-light radius (Jeft) or concentration 
(measured as R90/R50; right). The figure also shows that the correction should depend to first order on R e , whereas concentration is a second- 
order correction. The total sample of 8,00 galaxies is sp lit about the median value of R90/R50, showing the most concentrated galaxies in 
black. Our results are consistent wifh lGraham et alj d2005h . as concentration correlates with Sersic index. However our regime of observations 
is different from their work, as our galaxies have sizes significantly closer to the PSF. 

The simulations shown in figure [AT] motivate a simple correction dominated by R e with concentration modelled as a linear change to 
the zero point of the correction. The corrections - shown in the leftmost panels of figure [ATI for the median value of Rgo/Rso - are: 



with the radii given in pixels. These two expressions are used to make a correction on apparent magnitude and half-light radius. In order to 
assess the accuracy after the corrections we generate four samples of 910 simulated galaxies with the same apparent magnitude and size as 
our observed sample. The Sersic index, ellipticity and position angle are chosen randomly in the same way as described above. Figure lA2l 
shows the histograms comparing recovered and original size, magnitude and apparent surface brightness after applying the corrections given 
by equations lAll and lA^l The lines are Gaussian fits to the data, with the standard deviation shown in the top left corner of each panel. These 
are our official uncertainties, namely 9% in R e , 0.05 mag in the TOTAL iF775w magnitude and 0.16 mag/arcsec 2 in iF775w apparent surface 
brightness. 



APPENDIX B: RECOVERING STELLAR MASSES 

The stellar masses presented in this paper are extracted from the total apparent ip775w magnitudes. One could argue that NIR photometry 
might be a more reliable tracer of the stellar mass. To confirm the validity of iF775W photometry over the redshift range probed by our 
sample, we show in figure lBTl the predicted mass for a if775\v = 23 (AB) galaxy (top) or a K s = 22 (AB) galaxy (bottom). The result is 
shown according to photometric type, with the left (right) panels showing red (blue) galaxies. The models are identical to those described 
in §§2.3. The dark grey shading gives the expected stellar mass as a function of redshift, and - to guide the eye - we include a vertical line 
at z=0.7 (roughly the median redshift of the sample). As expected, for a later (i.e. bluer) photometric type, the stellar mass is smaller. One 
should emphasize here that within our volume-limited subsamples SI and S2, the fraction of galaxies corresponding to the panels on the 
right amount only to 4 and 11%, respectively. 

The horizonal light grey shading shows a 0.3 dex uncertainty in stellar mass for a z=0.7 galaxy. Notice that in this redshift range iF775w 
is as accurate a tracer of stellar mass as K 3 , notwithstanding the fact that the high resolution of the ACS images allows us to determine more 
accurately the TOTAL apparent magnitude of the galaxies. 

As a second check, we retrieved the publicly available ISAAC/VLT if s images of the GOODS/CDFS field (Retzlaff et al., in prepara- 
tion). The images were analysed with SExtractor dBertin& Arnoutsll 19961) using the provided zero point and weight maps. The catalogue of 
detected sources was cross-correlated with our sample, and only those galaxies with a SExtractor f lag=0 were selected for the comparison. 
We used MAG_BEST and further applied the same correction to obtain total magnitudes as in appendix A. Figure lB2l shows the comparison, 
along with the equivalent Gaussian distribution, which has a mean of (log M a (K s ) — log M s (i)) = —0.04 dex and standard deviation 
0.13 dex. The mean and standard deviation are reduced to —0.02 and 0.11 dex, respectively when the sample is restricted to zf775w < 23 
galaxies. We conclude that our i-band extracted stellar masses are accurate to within the quoted 0.3 dex uncertainty (including population 
synthesis effects). 
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Figure Al. Recovered ip77sw (top), and half-light radius (bottom) of simulated galaxies, lcr error bars shown. Each symbol corresponds to a choice of 
R e . The sample of 8,000 simulated galaxies is separated with respect to R90/R50, such that the grey (black) dots represent the subsample below (above) the 
median. The highest discrepancy in the recovery of the half-light radius occurs for values close to the FWHM of the PSF (~ 3 pixels). The simulations are 
shown with respect to Re (left), concentration middle - measured as R90/R50 - and apparent magnitude (right). The lines in the leftmost panels give the 
corrections used for flux and half-light radius (equations IA1 l and lA2t at the median value of concentration. Notice that the panels on the right only show the 
subsamples with the largest and smallest value of R e to avoid crowding the figure. 
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Figure A2. Recovery of four simulated samples of 910 galaxies with the same distribution in apparent magnitude and size as the original sample. The 
histograms compare the original and the recovered size (left); magnitude (middle) and apparent surface brightness along with a Gaussian that characterizes the 
uncertainty for the complete sample. 
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Figure Bl. Comparing the accuracy of stellar mass estimates using iF775w (top) vs a more "traditional" method with NIR magnitudes (bottom). A fiducial 
apparent magnitude is assumed as labelled. The dark grey region gives the stellar mass expected as a function of redshift for photometrically "red" {left) 
and "blue" galaxies. To guide the eye, we assume a z=0.7 galaxy (vertical line). The light gray horizontal shading represents a 0.3 dex uncertainty about the 
predicted value. Notice both ip775w an d Ks give similar accuracy over the redshift probed by our sample. 
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Figure B2. Comparison between stellar masses from this work - extracted from the corrected F775W magnitudes - and NIR masses using the K s magnitudes 
of the CDFS using the ISAAC images (Retzlaff et al., in preparation). Identical flux corrections are applied to the NIR magnitudes as those obtained for F775W 
(see appendix A). The grey histogram corresponds to a subsample of the brighter («F775w < 23) galaxies, which give an offset of = —0.02 and a standard 
deviation of a = 0.11 dex. 



